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Abstract 

In this paper we study the effects of s —* d + 2gluon on Kl — > 77 in the 
Standard Model. We find that this interaction can induce new sizeable SU(3) 
and U(3) nonet breaking effects in Kl — rj, rf transitions and therefore in 
Kl — > 77 due to large matrix elements of { r q{r}')\a s G < f lv G^ u \Q) from QCD 
anomaly. These new effects play an important role in explaining the observed 
value. We also study the effects of this interaction on the contribution to 
Am KL „ Ks . 

PACS numbers: 11.30.Hv, 13.25.Es, 14.40.AQ 
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It is well known that contributions from intermediate hadronic state effect play an im- 
portant role in many low energy processes. Some of the notable examples are K L — > 77 [1,2] 
and ArriK = ttik l — rriK s [3-7]. For Kl — > 77, the direct contribution due to quark level 
s — > c/77 alone accounts for only a small portion of the amplitude measured experimentally 
[2,7]. For Am K , the direct contribution due to AS = 2 four quark operator is again only 
a fraction of the experimental value depending on the value of the bag factor Bk [3,4]. A 
simple method to estimate the contributions from intermediate hadronic states is the pole 
dominance approximation in which one assumes that a few low lying resonances saturate 
the contribution. The commonly identified resonances in the above two cases are it , r\ 
and rj'. Combining with U(3) flavor symmetry, the K L — > 77 amplitude can be estimated 
[2,6]. If U(3) nonet is a good symmetry, the calculations are straightforward. However, 
not only nonet but also SU(3) are known to be broken, there are large uncertainties in 
these calculations. One should also study that if there are some new contributions in the 
Standard Model (SM) which have not been examined so far. In this paper we show that 
indeed there is a new contribution to to K L — > 77 and Am K . This new contribution comes 
from s — > d + 2gluon induced K — 77(77') transition, and the intermediate 77(77') subsequently 
decay into 77 or change to another neutral kaon through the usual As = 1 interaction. We 
find that, because of the large QCD anomaly hadronic matrix element (77(77') |a s G^G^|0) 
(G% v = e^ uaf3 G a a() ) , the new contributions are sizeable and also induce new sizeable SU(3) 
and U(3) breaking effects. 

The decay amplitude A^ ir of the direct contribution to — > 77 from quark level 
interaction s — ■> C/J77 in the SM has been studied before [1,2,7]. Here we improve the 
calculations by including QCD corrections which also serve to set up our notations. In the 
SM, s — > c/77 can be generated at one loop level by exchanging a W boson and quarks with 
two photons emitted from particles in the loop and particles in the external legs. The QCD 
corrected effective Hamiltonian for s — ■> c/77 is given by 

H eff (s - rf 7 7) = M]l + Ml\ (1) 
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where M]^ is the irreducible contribution with the two photons emitted from particles in the 
loop. M^ 7 is the reducible contribution with at least one photon emitted from an external 
s or d quark. 

The irreducible contribution M]^ is given by [2,7,8] 

= _ ^^G F n _1 ^ v * dVisF{X:Xi) d Y LR ^se^ih), (2) 

971 2f£ l ■ k 1 i=u,c,t 

Here e^{k) is the photon polarization vector with momentum k, L(R) = [1 — (+)7 5 ]/2, 
N = 3 is the number of colors, a 2 = C\ + c 2 /N, x = 2k\ • kijm^^ Xi = mf/m^r, and 

k 2^v P a\ k \ k \ + ^i ' k 2 e^ vp(T {k 2 - h) a ■ The function F(x,Xi) is given 

by 

F( X , Xi ) = (3) 

x Jo y 
The reducible contribution M^ 7 is given by [7,8] 

6tt j= ^ t p d -ki p s - k 2 p d -ki p s - ki 

+ (h -> fc 2 , fc 2 -> |/ -> i/, i/ -> A*)](m d L + m s R)se* ll {k 1 )e* v {k 2 ). (4) 

In the above q are the Wilson coefficients defined in the following AS 1 = — 1 effective 
Hamiltonian [9] 

AC 

H eff (AS = -1) = -^TOdd + c 2 2 ) -E E ^(40*)], (5) 

fc i=u,c,t 

where the summation over k is on all possible operators, four quark operators, quark-photon 
and quark-gluon operators, which are defined in Ref. [9]. The operators directly relevant to 
our calculations to the leading order are 

O x = q'y ll Lqd'fLs, 2 = d^Lqq^Ls, 

g _ 

07 7 = j^da^F^trridL + m s R)s, 

8G = -^da, u T a G^{m d L + m s R)s, (6) 

where G%" and F^ u are the gluon and photon field strengths. Here we have also written 
down the operator 0$g which is needed for the study of s — > dgg. 



To obtain the amplitude A dir for Kl — > 77 from the effective Hamiltonian H e ff(s — > 
c/77), one needs to bind the d and s quarks to form a kaon which involves long distance 
non-perturbative QCD effects. This effect cannot be calculated at present and is usually 
parameterized by a decay constant fx as, (0\d'y^'y 5 s\K ) = —if K P^- with f K determined 
from data. We have, 

A dir (K° -> 77) = (ry\H e ff(o 7j )\K°) 

= 2y/ ^ Gp f K [i(Na 2 V: d V us + 3^V td V ts )F, u F^ + 3^V td V ts F, u F^}, (7) 

where if e //(07 7 ) indicates the term proportional to Oj y in the effective Hamiltonian of eq. 
(5). F, v = (1/2) W F^. 

In obtaining the above result, we have used the fact that F(x,x Ci t) ~ —1/2 (large x c>t /x, 
and F(x,x u ) ~ (small x u /x). We also neglected small contributions from Cy 7 c which are 
proportional to x UjC [10], but have kept (* 7 which is —0.3 in the SM. 

The parameter £ is an average value of the quantity, k = —{m 2 K /l§){l/p d -ki — l/p s ■ k 2 + 
l/Pd'^2 - l/Ps - ^)- H one assumes that the 0? and s quarks share equally the kaon momentum, 
then £ = 1 [2]. We have also estimated £ by calculating the quantity < 0|/cd(l + 75)31^° > 
using perturbative QCD method and appropriate distribution amplitude of quarks in the 
kaon [11]. This approach also obtains a value of order one for £. One should be aware that 
the applicability of pQCD may not be a good one here. However, we find that contribution 
related to £ is not important as long as £ is of order one. That is, the precise value of £ is 
not important here and we will use £ to be one in our later discussions. 

To estimate the irreducible contribution, one needs to know the quantity a 2 = c\ + c 2 /N . 
Without QCD corrections, c± — and c 2 = 1. This gives a a 2 = 1/3. With QCD corrections 
the value for a 2 will be altered. The leading and next-leading order corrections to q have 
been calculated [9]. The values of q depend on the renormalization scale ji. Since one 
does not know precisely where is the matching scale //, this causes uncertainty in a 2 . For 
example at the leading order, a 2 = —0.27 at // ~ 1 GeV, while at /i = 1.3 GeV, a 2 = —0.17 
with Ajg-g = 325 MeV. At the next leading order the dependences on /i for each of the C\ 



and c 2 are reduced, but leaves a 2 still sensitive to fi. For example, in the NDR scheme, 
for = 325 MeV, a 2 is -0.08 and -0.1 at /i — 1.0 GeV and fi = 1.3 GeV, respectively. 
Allowing the QCD parameter A^r to vary within the allowed range 215 ~ 435 MeV, a 2 
can vary in the range —0.1 ~ —0.35 depending whether NDR or HV scheme is used [9]. 
That is, the value of a 2 is not well determined even from the next leading order perturbative 
calculations. When all effects, perturbative and non-perturbative, are correctly treated, the 
final physical observables will not depend on the renormalization scale ji. Unfortunately, 
such a calculation is not possible at present. The parameter a 2 behaves similarly to the 
one in hadronic B and D decays. In both D and B decays, the parameter a 2 determined 
from data (|a 2 | ~ (0.2 ~ 0.5)) is very different from factorization value by inserting ci j2 
at relevant scale in the expression for a 2 [12]. One would expect similar thing happens in 
kaon decays although the details may be different. To take into account uncertainties in 
theoretical calculations of a 2 , we will treat it as a free parameter and allow it to vary in the 
range of —0.5 ~ 0.5. One can also turn the argument around to obtain information about 
a 2 from K L — > 77 data. 

For £ of order one, and a 2 in the range of —0.5 ~ 0.5, we find that the dominant direct 
K° — > 77 amplitude is from the irreducible contribution. We have 

A dir (K L -> 77) = iAu^F^F^, 

A dir = ^^f K Na 2 Re(V: d V us ). (8) 

Using V ud = 0.9735 and V us = 0.2196 and f K = 1.27 f n [13], we obtain, 

A dir = 2.54 x lO^aaMeV- 1 . (9) 

For [02] = 0.5, it is only about 35% of the experimental value of 3.5 x 10~ 12 MeV -1 [13]. 
Without QCD corrections a 2 = 1/3, A dir is about 24% of the total amplitude. There must 
be some other contributions to this process. These effects may come from contributions 
with intermediate hadronic states or even contributions from new physics beyond the SM. 
If one has a good understanding of all SM contributions, one can have a detailed study of 
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new physics beyond the SM. It is probably too early to say that new physics is needed here 
due to large uncertainties in possible hadronic intermediate contributions. Therefore we will 
work within the SM and see how contribution from hadronic intermediate states can affect 
the results. 

Several analyses have been carried out using pole model with 7r°, t? and rj' poles to 
calculate the hadronic intermediate contribution. In this model, the amplitude Ahad from 
exchange of intermediate hadronic states is given by [6] 

(tt°\H w \K l ) 



A had = A(tt ^77)- 



m K ~ m l 



X 



to 



ml i(r]^77) (l + S , , 2^2 



1 + — £ l^r 1 , — —^cos6 + -^psinO 



m K — TO; 



2 i (7r o ^ 77 ) ^ ^ 73 

m\ - m\ A(i ^77) fl + S . „ 2^2 ^ 
H y y^T ~~r^ sm 9 -TP cos 9 



TO 



K 



m 2 v , A(n° -> 77) \ V3 



(10) 



where 9 is the r/ — rf mixing angle, 5 is the SU(3) breaking parameter [6]. The parameter p 
parameterizes U(3) nonet breaking effect and is defined as 

[3 <m\H w \K° > 
P V8< n°\H w \K° >• 1 J 

In the nonet limit p = 1. Chiral Lagrangian analysis gives (n \H w \K L ) = 1.4 x 10~ 7 m 2 K [6]. 
Using experimental values for 7r°, 77, 77' — > 77, Ay ia d can be estimated. 

The above contributions can be viewed as obtained by Ahad = J2i(ll\i){i\H e ff(AS = 
— 1)\K L ) with % = 7T°, 77, 77' in the pole model approximation. Here H e ff(AS = —1) is the 
full AS = — 1 effective Lagrangian with O^^ term removed since it has been counted as the 
contribution to A dir . Therefore A dir and A had are contributions from different sources. In 
previous calculations the contributions for A had from O s g were not considered [2,6]. We now 
study in detail the effect of this interaction on Kl — > 77. 

At the quark-gluon level, Osg induces s — * d + gg. To obtain Ah ad , one needs to estimate 
the contribution from s — > d+gg to K° — 77, rj' through gg — > 77, 77'. The effective Hamiltonian 
Mirr for s — > d(?<7, with color singlet ds bi-spinor product, can be obtained by some simple 
replacements from M] RR . To obtain Mf R R one first replaces the photon polarization vectors 
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e M (/ci)e^(/c 2 ) by the gluon polarization vector e^(/ci)e„(/c 2 ) with the color index a summed 
over. Then one replaces a em by a s (9/4)/ (27V) and a em c[ 2 by a s c\ iG /{2N) for Mj 9 R and Mff , 
respectively [7]. The factor 1/(2N) comes from picking up the color singlet part. 

Similar to the procedure in obtaining the amplitude A dir for K L — > 77, one can obtain 
the amplitude for K L — > gg. We find that with £ of order one, a 2 in the range of —0.5 ~ 0.5 
and Cg G ~ —0.15 as given in the SM, the irreducible contribution, again, dominates the 
amplitude. We have 

A(K L -> gg) = ^ 2 -^f K Na 2 Re(V: d V us Y\G%G^. (12) 

The above interaction can induce large Kl — 77, 77' transitions and therefore contribution 
to K L — ► 77, because QCD can induce large matrix elements for {r]{ri')\a s G a G^ v \Qi) . 

QCD anomaly implies that the divergence of the singlet current, a 1 = ^7^75^ + ^7^756?+ 
57^755, is not zero in the limit of zero quark masses, and is given by 

(77(77') |d M a*|0) = (r)(r)')\2i(m u u7 5 u + m d d^ 5 d + m s s-f 5 s)\0) 

- (^OI^^Gf |0>. (13) 

While for the octet current, = vry^^u + d^^/y^d — 2s7 M 7 5 s, one obtains [14] 

(77(77') |9"aJ|0) = (77(77')|2i(m u iZ75M + m^d - 2m s s 75 s)|0). (14) 

Since m u rf are much smaller than m s , one can neglect terms proportional to m U}d . One then 
obtains 

W(P)\^GI U G^\0) = ^(v^cos^/gsin^ 2 , 

(V(P)\^G%G^\0) = ^(-v^sin^/gcos^ 2 , (15) 

where /i j8 are the singlet and octet pseudo-scalar decay constants. 

If there is no 77 — 77' mixing and all quark masses are equal, the gg state being a flavor 
singlet can only have transition to 771. Because the 77 — 77' mixing and the different quark 
masses, both U(3) nonet and SU(3) symmetries are broken. The Kl — > 77,77' transitions 
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induced by s — > dgg will induce nonet and SU(3) breaking in the total amplitude A total . 
Normalizing the signs of each contributions to theoretical calculations, we finally obtain 



" m\ - ml A(r) -> 77) + 5 + 2^, 0(A . ^ 

x H f LLL cos6» + -^(p + r 9S )sm6» 



mt- — m: 



A' 



i(vr° -> 77) V 



~ ^ ^ " T 'M iTU ; - sin g - ^(p + r") cos g 



2 



m 2 A(rf ^ jj) (1 + 5 + 539 . 2>/2, 



(16) 



where 5 9S and r" are the SU(3) and nonet breaking induced by the s — > dgg interaction. 
They are given by 

xgg f _2 G FRe(V* d V us ) „ 

6 = -V2f K f 8 m K 02, 

= _ A^. (17) 

We find 

593 = o.96 ^a 2 , r" = -0.48^a 2 . (18) 
Ik Jk 

We see that the corrections can be sizeable and can not be neglected. 

We now provide some details for numerical calculations. There are several parameters 
involved in Ahad, the mixing angle 9, the decay constants fi t $, the SU(3) and U(3) nonet 
breaking parameters 5 and p, and the parameter a 2 . Chiral perturbation calculations and 
fitting data not involving K L — > 77 have obtained 9 rs —20°, 5 0.17, / 8 p=* 1.28/^ 
and /1 ~ l.lO/n- [15]. We will use these values for these parameters in the calculation of 
Kl — > 77. There is not a reliable estimate for the parameter p. Since we are interested to 
see how the new s — > dgg interaction induces U(3) nonet breaking effect, we will take p = 1 
and attribute nonet breaking solely to r" . As have been discussed s — > dgg also induce 
SU(3) breaking effect. This effect was not included in other fittings. We therefore should 
include this new SU(3) breaking effect also. 

Without the s — > dgg effect, we find that the amplitude A total is equal to 5.5(l + 0.46a 2 ) x 
10~ 12 MeV -1 which is considerably larger than the experimental value 3.5 x 10~ 12 MeV -1 
[13] for |a 2 | < 0.5. With the new effect, we find 
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total 



5.5(1 + 2.14a 2 ) x 10" 12 MeV 



-i 



(19) 



To reproduce the central experimental value, a 2 is required to be —0.17 which is a reasonable 
value to have. 

The detailed numerical results depend on several parameters. Even with other parame- 
ters fixed, one can introduce also a phase to a 2 . To fit the K L — > 77 data, the values for the 
magnitude and phase of a 2 can vary. We, however, would like to emphasize that the new 
effect discussed can play an important role in K L — > 77 independent of the details. 

The new contributions for Kl — 77(77') transitions also induce new hadronic intermediate 
state effect to the K L and K s mass difference parameter Re(M 12 ) in the pole dominance 
approximation. We find [6] 



non-negligible portion of the experimental value of 3.5 x 10~ 12 MeV. With the new effects 
and a 2 = —0.17 as determined from Kl — > 77, the contribution to Aijik is —0.9 x 10~ 12 
MeV, and again it can not be neglected. The new effect in K L — > 7r°, rj, rj' transitions can 
have sizeable contribution to Atuk- 

The s — > dgg process can also induce i^L-glueball mixing, which would also affect Kl — > 
77 and Ams-L, as pointed out in Ref. [7] where a light glueball mass 1.4 GeV was used. 
Recent lattice calculations indicate that the pseudo-scalar glueball mass is about 2.3 GeV 
[16]. With such a large mass the glueball-?] (r/) mixing contribution should be small and 
therefore the effects are smaller than effects discussed earlier. 

In conclusion we have evaluated additional contributions to Kl — > viv') transitions from 
s — > dgg in the Standard Model. These transitions induce new sizeable SU(3) and U(3) 
breaking effects and have significant effects on contributions to Kl — > 77 and Aijik- 



2m K Re(M 12 ) 



\{AHw\K») 




(20) 



Without the new effects, the above would lead to Am K = —0.5 x 10 12 MeV which is a 
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